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1 Introduction 



Nuclear shadowing is a well-established phenomenon, which attracts attention of both ex- 
perimentalists and theoreticians. It was found (see [U [5] and references therein) that the 
inclusive nuclear structure function is smaller in nuclei than in a free nucleon at small values 
of the Bjorken variable x < 0.01. The nature of shadowing, emerging e.g. in deep inelastic 
scattering (DIS), can be well understood in terms of multiparticle scattering in the target rest 
frame. The incoming photon is represented as a superposition of gluons, quarks, anti-quarks 
and their bound states. At high bombarding energies the photon converts into qq-pair long 
before the target, and its hadronic component interacts coherently with several nucleons of 
the target nucleus. This process leads to absorption and, therefore, to nucleon shadowing 
(for a review see e.g. [3]). Among the important consequences of the phenomenon is, for 
instance, severe reduction of particle multiplicity in heavy-ion collisions at LHC energies 
(-y/s = 5.5 TeV), since multiple scattering is connected to diffraction [HElll]. The effect 
can be further decomposed onto the quark shadowing and the gluon shadowing; the latter 
provides largest uncertainties in the theory and is a subject of our present study. 

In recent years a lot of interest has been generated about the possibility of parton sat- 
uration in the nuclear wave function at the smallest x accessible at HERA, and there is an 
ongoing discussion if the same effects can be observed for hadron-nucleus and nucleus-nucleus 
collisions at RHIC. In this paper, we will focus mostly on the low- a; effects mentioned above. 
Understanding the so-called cold nuclear effects in hadron-nucleus collisions serves also as a 
baseline for the correct treatment of possible final state effects in nucleus-nucleus collisions, 
e.g. high-pT particle suppression and heavy-flavor production at RHIC. 

Our starting point is noticing that a significant change in the underlying dynamics of 
a hadron-nucleus collision takes place with growing energy of the incoming particles. At 
low energies, the total cross section is well described within the probabilistic Glauber model 
[7] , which only takes into account elastic rescatterings of the incident hadron on the various 
nucleons of the target nucleus. Elastic scattering is described by Pomeron exchange. At 
higher energies, E > Ecrit ~ niNlJiRA (a* is a characteristic hadronic scale, /i ~ IGeV, and 
Ra is the radius of the nucleus) corresponding to a coherence length 



the typical hadronic fluctuation length can become of the order of, or even bigger than, 
the nuclear radius and there will be coherent interaction of constituents of the hadron with 
several nucleons of the nucleus. The sum of all diagrams was calculated by Gribov [SI [9], 
which corrected the Glauber series by taking into account the diffractive intermediate states 
in the sum over subsequent rescatterings. The space-time picture analogy to the Glauber 
series is nevertheless lost, as the interactions with different nucleons of the nucleus occurs 
nearly simultaneous in time. The phenomenon of coherent multiple scattering is referred to 
as shadowing corrections. 

An additional effect which comes into play at high energies, is the possibility of interac- 
tions between soft partons of the different nucleons in the nucleus. In the Glauber-Gribov 
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model this corresponds to interactions between Pomerons. These diagrams are called en- 
hanced diagrams ^U\, and can also be understood as interactions between strings formed in 
the collision. Actually, the necessity to include such diagrams at high energies can be related 
to unitarization of the total cross section. There is a connection between these effects and 
saturation effects already mentioned earlier. 

The Glauber-Gribov model is described and a unitarity-conserving procedure for finding 
the total cross section of 7*-nucleus (7*^) interaction, which corresponds to summing up 
Pomeron fan-diagrams, is presented in Sec. [H Further we will concentrate on new and 
interesting data on gluon diffractive distribution function, which we will describe in Sec. [3l 
The results for gluon shadowing are presented in Sec. 21 and our conclusions are drawn in 
SecE 

2 The Model 

We consider the nucleus as a set of nucleons, in the spirit of the Glauber model. The elastic 
J* A scattering amplitude can then be written as the sum of diagrams shown in Fig. [T] 
i.e. as multiple 7*-nucleon (7*iV) scattering diagrams with Pomeron exchange [HEI]. The 
contribution from 1,2... scatterings 

should be summed up to obtain the total cross section. In Eq. Q , the first term simply equals 
to the Glauber elastic contribution and subsequent terms describe multiple interactions of 
the incoming probe with the nucleons in the target nucleus. 

The multiparticle content of subsequent diagrams in Fig. [T] is given by AGK cutting 
rules [12], where the intermediate states are on-shell. The cut contribution of the double 
rescattering diagram can be expressed in terms of diffractive deep inelastic scattering (DDIS). 
The usual variables for DDIS: Q^,x,M'^ and t, or xjp, are shown in Fig. [51 The variable 
(3 = Qi^j^i — x/xp plays the same role for the Pomeron as the Bjorken variable, x, for the 
nucleon. We assume that the amplitude of the process is purely imaginary; this is justified 
for a value of the Pomeron intercept close to unity [13 . The contribution from the second 
term in Eq. ^ to the total 7* A cross section is given by [3] 

4'^^ = -A^A{A - 1) X 

t=0 

where = f^^dzpA{b,z) is the nuclear normalized density profile, J d^bTA{b) ~ 1. 

The form factor Fa is given by 

FAitrmn) = j A% M^/^t~b) T A{b) , (4) 

where tmm = —"mf^Xp, and Jq{x) denotes the Bessel function of the first kind. Strictly 
speaking, Eq. ^ is valid for nuclear densities which depend separately on b and z, however 
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we have checked that calculations with an exact expression lead to negligible corrections. 
Note that since Eq. ^ is obtained under very general assumption, i. e. analyticity and 
unitarity, it can be applied for arbitrary values of provided x is very small. We have 
assumed R\ 3> R%, so that the t-dependence of the '-f*N cross section has been neglected. 
For a deutcron, the double rescattering contribution has the following form 

where Foit) = exp(ai), with a = 40 GeV"^ g]. 

In Eqs. ^ and (O, M^j„ corresponds to the minimal mass of the diffractively produced 
hadronic system, M^j„ = 4m^ = 0.08 GeV^, and M^^^ is chosen according to the condition: 
xip < x^"^. The choice of x"^""^ is governed by the fact that the model is only valid for 
Xjp <C 1, i. e. a large rapidity gap is required in the experimental data. We use the standard 
choice for Xp""^ = 0.1 [13]. It is convenient as it guarantees the disappearance of nuclear 
shadowing at a; ~ 0.1 as in experimental data. Coherence effects are taken into account 
through FA{tmin) in Eq. (H]), which is equal to 1 at x ^ and decreases with increasing x 
due to the loss of coherence for x > Xcrit ^ (it^nRa) ^- In the numerical calculations, a 
3-parameter Woods-Saxon nuclear density profile with parameters from |15] has been used. 

Higher order rescatterings in Eq. ([2]) are model dependent. We use the Schwimmer 
unitarization '16] for the total cross section which is obtained from a summation of fan- 
diagrams with triple-Pomeron interactions. It was checked in g] that it gives results very 
close to other reasonable models, such as the quasi-eikonal model. The total cross section is 
then 



^Sch _ „ 



/• 2. ATA{b) 

J 1 + {A-l)f{x,Q^)TA{b)' 



where 



f{x,Q') = ^ / dM' 



dM2 At 



Equation ^ does not take into account the shadowing effects of valence quarks nor anti- 
shadowing effects, which may play an important role for x>0.1 [l]. 

Nuclear shadowing is studied in terms of the ratios of cross sections per nucleon for 
different nuclei, defined as 

R{AIB) = f ^ , (8) 

as a function of x, which can in turn be expressed via structure functions of the different 
nuclei. The simplest case is B=:N, then 

In this framework shadowing can also be calculated at a fixed value of the impact parameter 
b for given values of {x, Q^} 

i?^^'' (A/N) {b) = I . (10) 

1 + (A - 1)/(x,Q2)Ta(6) 
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Thus the total 7* A cross section can be calculated within the Glauber-Gribov model in a 
parameter- free way provided the total 7* TV cross section and the differential cross section for 
diffractive production are known. 



3 Diffractive parton densities from HERA 

The cross sections of the inclusive and diffractive processes are expressed through the nucleon 
structure functions, which are in turn associated with distribution functions of partons in 
the nucleon and in the Pomeron. In DIS the structure function of a nucleon, F2{x, Q^), is 
related to the total cross section of 7*iV interaction through 

0'7*JV = — — F2[X,(4 ) , 

valid at small x. Similar to the inclusive DIS case, a factorization theorem has been proven 
in perturbative QCD to hold for diffractive structure functions [T7]- The relation between 
the diffractive cross section and the diffractive structure function is given by 



dMMt 



xpF^^lQ ,xip,P) 



02 (g2 + M^) 



where the usual factorization has been assumed: 



da^.r^ix,Q^,M^,t) 
dM2 dt 



da^,j,{x,Q\M^ 



dAP dt 



J t=o 

A further assumption about the lower part of Fig. [21 the so-called Regge factorization [IS] , 
allows us to write the diffractive structure function as 

F^^{xjp,Q^P) - fjp{xjp)F{p,Q^) . (11) 

The first factor is referred to as the (t-integrated) Pomeron flux (we will not take into account 
the sub- leading contributions from other reggeon trajectories), which is in turn defined as 

fjp{xp) = Ap __^-^di , 

•^icut Xjp 

where Ap is a constant that fixes the normalization of the flux, and tcut — — IGeV^. As 
usual, we assume a linear Pomeron trajectory, ap (t) = ajp(O) -I- a'jpt. The second factor in 
Eq. pTjl . F{P,Q^), is the Pomeron structure function. For the sake of simplicity, Eq. ([7]) 
can now be reduced to 

fix,Q') = 4^ r dxpB{xp)^^^^^l^Flit^,^.), (12) 
Jx F2[x,Q^) 

where B[xp) = B^ + a^ln^. The ratio of structure fmictions in the integrand can be 
understood as the density of partons in the Pomeron compared to the density of partons in 
the nucleon. 

The determination of the distribution of partons in the nucleon and Pomeron cannot 
be carried out in pQCD, since it depends on soft processes such as conflnement. The in- 
clusive distributions are measured in DIS experiments to high accuracy and predictions of 
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the DGLAP evolution are in good agreement with experiment for a broad range of x and 
values. Yet, the experimental status of the diffractive distribution functions was, until 
recently, uncertain. The biggest uncertainty was related to the gluon distribution, because 
it is not measured directly in the experiment. Some authors therefore assumed that the 
shadowing from quarks and gluons are of equal strength [19l [201 EI] , while others predicted 
that the gluon shadowing dominates [HI [531 [H] • 

The results of new high-precision measurements of the diffractive parton distribution 
functions (DPDFs) presented by the HI Collaboration |2S1I2S], shed new light on the role of 
gluon shadowing at intermediate Q^. The quark singlet and gluon DPDFs were fitted by a 
simple function 



where Ai, Bi and Ci are fitting parameters, at a given Qq ~ 2 GeV^. In general, to ensure 
that the r.h.s. of Eq. (fTij) always disappears as (3 1, it should be multiplied by exp ( ) . 



QCD evolution to arbitrary was performed by the HI group. For the quark singlet 
distribution, all three parameters {Aq, Bq,Cq} were used in the fit, while the gluon density 
is found to be insensitive to the Bg parameter, which is therefore set to zero. This fit is 
referred to as the 'FIT A'. Additionally, because of DGLAP evolution, data at /3 > 0.3 cannot 
constrain the gluon density because of the large quark contribution to the evolution. This 
lack of sensitivity is confirmed by repeating the fit with the parameter Cg set to zero, and is 
further referred to as the 'FIT B'. The values of the corresponding Pomeron parameters are 
listed in Tab.[T] The fitted quark singlet distributions of 'FIT A' and 'FIT B' arc consistent 
with each other. 

Inclusion of diffractive di-jet production in the analysis provides a big improvement in 
the precision of the gluon density measurement \Tll , which can be further constrained by the 
diffractive production of charm in DIS [28j. A combined fit to DDIS data and diffractive 
di-jets [57j results in a curve similar to fit B (described above) yet with a slightly smaller 
gluon density at (3 > 0.5. For completeness, we also compare the new results with the old 
HI parameterization presented in 2002 (corresponding Pomeron parameters are given in 
Tabic [U. 

The gluon distribution is almost a factor of 10 bigger than the quark distribution at the 
same for a broad region of /3 |25l [26] . Shadowing from quarks was calculated in [11] using 
the old HI parameterizations, and we have checked that the new HI data are consistent with 
their calculations. Furthermore, in the relevant kinematical range for hadron-nucleus and 
nucleus-nucleus collisions at RHIC and LHC, the gluon density dominates. In what follows, 
we will therefore consider structure functions of gluons in nuclei. The first term in Eq. ([2]) 
is then proportional to AGn{x,Q'^), while the second rescattering term is correspondingly 



(1 - fif 
Agr""' (1 - /3) 



(13) 
(14) 




equal to -A{A - 1)Gn{x, Q^) J d^TA{hf f{x, Q^), where 




/3g^(/3,Q^) 




(15) 
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Thus the gluon shadowing factor, Rg, is calculated according to Eq. ^ with f{x, Q"^) given 
by Eq. IjlSp. The gluon distribution of the nucleon, Gn{x, Q^) = xg{x, Q^), was taken from 
CTEQ6M parameterization [3D]. 

The model for diffractive production by virtual photon described above is applicable 
for intermediate > 1 — 2 GeV^ [251 ESI ■ It is also important to keep in mind that the 
parameterization of HI data leads to a violation of unitarity for x — > and should be modified 
at very low x (probably already at a; '--^ 10^^ for low Q^) [31 . Therefore our predictions are 
reliable in the region x > 10^"' which is relevant for RHIC and most experiments at LHC, 
and should be taken with care for x < 10~*. 

4 Numerical results 

Gluon shadowing for deuteron and various heavy ions (Ca, Pd and Au) at — 6.5 GeV^ 
is presented in Fig. [3] We have made the calculations using two fits of the gluon diffractive 
distribution function, FIT A and FIT B. The gluon shadowing is very strong at small x, and 
disappearing at a; = x^"^. This is a consequence of the coherence effect in the form factor 
Eq. and the vanishing integration domain in Eq. Gluon shadowing is as low as 

0.35 for the Au/nucleon ratio. In Fig. d] shadowing in Pb is shown for different values of Q^. 
We have also included the result using the 2002 HI parameterization for — 6.5 GeV^. 
Although the difference for the experimentally measured gluon distribution function is quite 
substantial, the shadowing contribution is stable, changing maximally 30% for the lowest x. 

There is little change to be observed in the gluon shadowing ratio with increasing Q^. 
The QCD evolution of the main term and rescattering terms are effectively treated separately 
in our approach, and therefore the shadowing correction has a slow, logarithmic dependence 
on [35] . The dependence of shadowing was studied within a similar model in [3] and 
good agreement with existing data was found. 

A comparison of the results of our model with other predictions is presented in Fig. 
for Pb at = 5 GeV^. For all x < 10^^ our model predicts slightly stronger gluon 
shadowing compared to FGS [H] , while models based on global fits to existing data on nuclear 
modifications and DGLAP evolution [331|31] predict a modest gluon shadowing effect. The 
authors of [24j have made calculations in a similar framework as the presented model, using a 
quasi-eikonal summation of higher-order diagrams. In general, the quasi-eikonal model gives 
rise to a stronger shadowing effect than models including enhanced diagrams, i.e. Pomeron 
interactions, like the Schwimmer model, at x < 10^"^. This is not so clear from the figure as 
the authors of [21] put x^°-^ = 0.03 and include anti-shadowing for gluons. The quasi-eikonal 
model will, at asymptotical energies, lead to the grey disc limit. 

We have also calculated the impact parameter dependence of Glauber-Gribov gluon shad- 
owing given x = 10^^ and — 5 GeV^, and compared it to results from [21] in both "high" 
and "low" gluon shadowing mode. Our model predict shadowing in agreement with the 
strongest FGS prediction, which is 20% lower than in the "low" gluon shadowing mode 
for central impact parameters. The shape seems to be compatible between the two mod- 
els. We have also checked that agreement with the FGS model improves when updated HI 
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parameterizations are used |35j . 

5 Conclusions 

Calculations of gluon shadowing in heavy-ions and deuteron have been made within the 
Gribov- Glauber model, which is a parameter-free framework to calculate shadowing effects 
in hadron-nucleus collisions at high energy. Diffractive structure functions have been pa- 
rameterized using the latest NLO QCD fits from HERA experiments. Rather substantial 
shadowing effects are predicted for heavy nuclei at x < 10^^, which is slowly changing with 
increasing Q^. This will have implications for both light and heavy particle production 
in hadron-nucleus collisions at RHIC and LHC, where the relevant {a:,Q^} range is being 
probed. 

Comparison to other models predicting strong shadowing effects a.t x < 10^^ have been 
performed. Calculations of Glauber-Gribov shadowing using Schwimmer summation of fan 
diagrams differ from other approaches within the same framework mainly due to differences 
in the treatment of higher order rescatterings, which is important in the low- a: region, and 
due to differences in the QCD evolution in log Q^. The advantages of our method is the clear 
treatment of these effects within a well-established framework. 

The nuclear PDFs can be measured in ultra-peripheral collisions at both RHIC and 
LHC [35]. In the future, hopefully an electron-ion collider (eRHIC) would provide clear 
information about DIS on nuclei which would dramatically improve our understanding of 
high-energy nuclear effects such as shadowing. 
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of Figures 

The single and double scattering contribution to the total 7*A^ cross 

section 

DDIS kinematical variables in the infinite momentum frame 

Gluon shadowing for deuteron (dash-dotted curves) and for heavy 
ions: Ca (dotted curves), Pd (dashed curves), Au (solid curves). Black 

curves are for FIT A and grey curves are for FIT B 

Gluon shadowing for Pb for different virtualities, 

Comparison of the results of the Glauber-Gribov model with FGS 

model |24j . EKS [33] and HKM [Mj parameterizations 

Impact parameter dependence of gluon shadowing in the Glauber- 
Gribov model, compared to the predictions of FGS [23] 
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Figure 3: Gluon shadowing for deuteron (dash-dotted curves) and for heavy ions: Ca (dotted 
curves), Pd (dashed curves), Au (sohd curves). Black curves are for FIT A and grey curves are 
for FIT B. 
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Figure 4: Gluon shadowing for Pb for different virtualities, . 
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Figure 5: Comparison of the results of the Glauber- Gribov model with FGS model [23], EKS 
[33| and HKM [31] parameterizations. 
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Figure 6: Impact parameter dependence of gluon shadowing in the Glauber-Gribov model, 
compared to the predictions of FGS 
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List of Tables 

1 Pomeron parameters from HI Collaboration [25\ [26| [29] 
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ajp(O) [FIT A] 


1.118 


ajp(O) [FIT B] 


1.111 


a'jp 


0.06 GeV-2 


Bo 


5.5 GeV-2 


ajp(O) [2002] 


1.173 


a'jp [2002] 


0.26 GeV-^ 


-Bo [2002] 


4.6 GeV-^ 



Table 1: Pomeron parameters from HI Collaboration |25i 1261 [29 ] . 
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